Magneto-electric equivalence and emergent electrodynamics in bilayer graphene 
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It is a fundamental paradigm that the physical effects induced by electric fields are qualitatively 
different from those induced by magnetic fields. Here we show that electrons at a Dirac point in 
bilayer graphene experience an unusual type of electromagnetism where magnetic and electric fields 
are virtually equivalent: every coupling of an electron's degrees of freedom to a magnetic field is 
matched by an analogous coupling of the same degrees of freedom to an electric field. This counter- 
intuitive duality of matter-field interactions enables novel ways to create and manipulate spin and 
pseudo-spin polarizations in bilayer graphene via external fields and leads to the emergence of a 
valley-contrasting axion electrodynamics, where the traditional association of charges at rest with 
electric fields and charge currents with magnetic fields is reversed. 



I. INTRODUCTION 

The dynamics of charge carriers in crystalline solids 
normally resembles that exhibited by free electrons, only 
with numerical values of parameters such as mass and 
gyromagnetic ratio renormalized due to the influence of 
the crystal structure.^ In certain instances, changes in 
the properties of band electrons tu rn o ut to be dra- 
matic. Few-layer samples of graphite^' are one exam- 
ple, and a new class of so-called topological ma terials^ 
another. Here we focus on bilayer graphene^ (BLG), 
i.e., two single-layer sheets of graphite stacked as shown 
in Fig. [ija). Electrons in this material turn out to be 
pseudospin-carrying chiral fermions having a finite band 
mass but zero rest energy.^ Our study shows that the 
interaction of these exotic charge carriers with electro- 
magnetic fields is very unusual. Normally, electric fields 
couple to electric charges at rest and moving magnetic 
moments, whereas magnetic fields couple to moving elec- 
tric charges and magnetic moments at rest. In contrast, 
every coupling of a bilayer electron's degree of freedom 
to an electric field is matched by an analogous coupling 
to a magnetic field. We explore physical consequences 
of this magneto-electric equivalence, including anoma- 
lous polarizations of real spin and the layer-index-related 
pseudospin, as well as an anisotropic version of axion 
electrodynamics.^ 

The BLG band structure near the K point in the Bril- 
louin zone [see Fig. |ljb)] is described by the effective 
Hamiltonian 

— (fc_ er + + k + er_) , (1) 

where h is Planck's constant, k = (k x , k y ) is the electrons' 
wave vector measured from K, and the Pauli matrices 
&x.y,z are associated with the sublattice (or, equivalently, 



the layer-index) pseudospin degree of freedomP In our 
notation, <7o is the 2x2 unit matrix, a± = (a x ±ia y )/2, 
and k± = k x ± ik y . Numerical values fo r th e effective 
masses n%j and the speed V3 are well knownPSl Very close 
to the K point, the energy dispersion resulting from ([I]) 
mimics that of massless Dirac electrons, as is the case in 
single-layer graphene (SLG). However, as mi -C rri2, the 
dominant behavior of electrons in BLG is captured by 
the quadratic dispersion shown in Fig. [ljc) . 

External fields can have striking effects on the elec- 
tronic properties of charge carriers in solids. Previously, 
only the effects of electric fields £ z and magnetic fields 
B z directed perpendicular to the BLG sheet have been 




FIG. 1: Basic structural and electronic properties of bilayer 
graphene. (a) Honeycomb structure of a bilayer stack of 
graphene. Atoms in sublattice A (B) are marked in grey (red). 
A yz plane is marked in light grey, (b) Brillouin zone and its 
two inequivalent corner points K and K'. The remaining cor- 
ners are related with K or K' by reciprocal lattice vectors, 
(c) Dispersion E(k) near the K point. We have k = ft K. 



2 



TABLE I: Magneto-electric equivalence of matter-field interactions in BLG. In each row, the coupling of an electron's degrees 
of freedom to a magnetic field B (left) is matched by an analogous coupling to an electric field £ (right). In these pairs, B and 
£ couple identically to the real spin s and pseudospin a, while the coupling to the valley isospin r is opposite: For one field, 
the coupling has the same sign in both valleys (isospin to), for the other field it has opposite signs (isospin r z ). In each row, the 
interaction unique to BLG is marked by * , while the other one exists also in SLGP The interactions are visualized in Fig. [2] 



magnetic field B 
orbital Zeeman splitting (_L field) B z a z t z 

magnetic spin splitting (_L field) B z s z to 

magnetic spin splitting (|| field) (B x s x + B y s v )to 

spin-orbital Zeeman (_L field)* B z (s y a x to — s x a y t z ) 
spin-orbital Zeeman (|| field)* s z (B y a x tq — B x a y t z ) 
orbital Zeeman splitting (|| field)* (k x B y — k y B x )a z to 
trigonal Zeeman splitting (s y B y — s x B x )a x to 

+ (S X By + Sy B X )Oy T Z 



considered, giving the extended effective Hamiltonian 
K K (k ) ^,B z )=K K (k + eA) 

+ y £ z a z - I fi B B z a z + c z £ Z B Z a . (2) 

Accordingly, (i) a potential difference between the two 
lay ers ( equivalent to finite £ z ) opens up a pseudospin 
g ar j6iiQ| \ z £ z ^ (jj) q z induces a pseudospin Zeeman split- 
ting g/i B B z , and (iii) the simultaneous presence of 
fields £ z and B z leads to an (actu ally , valley-contrasting 
- see below) overall energy shiftPJK3c z £ Z B Z . In Eq. ([2]), 
g is the gyromagnetic ratio, /ib is the Bohr magneton, 
and A is the electromagnetic vector potential satisfy- 
ing (V x A) z — B z . The matter-field interactions (i)- 
(iii) generate sizeable effects for typical values of £ z and 
B z (see Appendix [C]) . In addition, the chiral nature of 
electrons in BLG leads to unconventional quantum-Hall 
physics 

II. MAGNETO-ELECTRIC EQUIVALENCE IN 
BLG: ORIGIN & PHYSICAL CONSEQUENCES 

Inspection of Eq. ([2]) reveals a surprising feature: dis- 
regarding constant prefactors, the electron's interaction 
with fields £ z and B z is symmetric with respect to the 
interchange of £ z and B z . Indeed, this observation is 
not an accident. It reflects the unusual property of BLG 
that its crystal symmetry does not distinguish between 
polar vectors such as the electric field £ and axial vec- 
tors such as the magnetic field B. Moreover, the familiar 
constraints due to time-reversal invariance are modified 
in BLG such that symmetry under time reversal likewise 
permits that £ an d B be come interchangeable. 

Using symmetry * 9 * 14 * 15 ! we have obtained the complete 
set of interactions of electrons near the K point in BLG 
with magnetic and electric fields B and £, strain e, and 
also spin s. See Appendix [A] for details on our method 
and full results. All interactions that exist in SLG are 
also present in BLG. In addition, there are new interac- 
tions in BLG that can be obtained from field-dependent 



electric field £ 

(1) £ z a z to inter-layer (pseudo-spin) gap* 

(2) £ z s z t z electric spin splitting (_L field)* 

(3) (£ x s x + £y Sy)T z electric spin splitting (|| field)* 

(4) £ z (s y a x t z — s x Oy to) Rashba spin splitting (_L field) 

(5) s z (£ y a x t z — £ x G y to) Rashba spin splitting (|| field) 

(6) (k x £ y — k y £ x )o z t z orbital Rashba splitting 

(7) (s y £ y — s x £ x )a x t z trigonal Rashba splitting* 
+ {s x £ y + s y £ x )a y to 



SLG interactions by replacing the magnetic (electric) 
field components with their electric (magnetic) counter- 
parts. Table [I] summarizes this magneto-electric equiva- 
lence for couplings that are linear in electric or magnetic 
field components. The interaction oc £ z a z coupling the 
pseudospin to a perpendicular electric field thus arises by 
magneto-electric equivalence from the pseudospin Zee- 
man term oc B z a z that exists in SLG. Moreover, we 
obtain a rather counter-intuitive purely electric-field- 
dependent spin splitting cx £ z s Zl which can be directly 
measured using electron spin resonance. Also interesting 
is the spin-orbital Zeeman splitting oc B z (s y a x — s x a y ), 
which is the magneto-electric equivalent of Rashba spin 
splitting in SLGP^ We note that the electrons interact 
not only with the fields B and £, but their quantum dy- 
namics depends also - via the usual minimal coupling 1 
- on the electromagnetic potentials A and $. The lat- 
ter contributions are not affected by the magneto-electric 
equivalence. 



A. Spin textures 

The interactions shown in Table [I] give rise to a diverse 
set of textures for the induced spin and pseudospin ori- 
entations in k space, which are illustrated in Fig. [2j The 
grey arrow in the lower left of each panel indicates the 
direction of the external field while a large arrow in the 
lower right indicates (if present) the spin polarization in- 
duced by the particular interaction when averaging over 
all occupied states in the k x k y plane. 

In addition to the layer-index pseudospin <x, electrons 
in BLG carry a valley-isospin r that distinguishes states 
near the two inequivalent K and K' = — K points in the 
Brillouin zone 3 [Fig. [ljb)]. The effective Hamiltonian for 
electrons in the K' valley can be obtained from that for 
the K valley by a reflection of the vectors k, s, £, and 
B at the yz plane, see Fig. [jJa)P Choosing the conven- 
tion that the B and ^-dependent terms have the same 
sign in the K valley (first column of Fig. [2| , the corre- 
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valley K : B, £ valley K' : B valley K' : £ 




FIG. 2: Visualization of real and pseudo-spin textures in BLG generated by field-dependent interactions shown in Table [I] 
[omitting the trivial interactions (2) and (3)]. For selected points in k space, a red (green) arrow indicates the expectation 
value of the real (pseudo-)spin vector of the corresponding negative-energy eigenstate obtained by diagonalizing the leading 
field-independent contribution to the BLG Hamiltonian together with the interaction associated with the same number in 
Table [i] The 1st (2nd, 3rd) column shows results for the K (K' with B, K' with £) valley. Note that our sign convention 
is such that B and £ have the same effect in the K valley (1st column). Grey arrows indicate the direction of the applied 
field. The first row denoted (0) shows the pseudo-spin orientation without any external fields. Big arrows in the lower right of 
the panels indicate the existence of a net spin polarization obtained by averaging over the (pseudo-)spin expectation values of 
occupied hole states for a negative chemical potential. 
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sponding term in the K' valley involving the axial vector 
B (second column) differs by an overall minus sign from 
the term involving the polar vector £ (third column of 
Fig. [2j see also Table [I]). For example, the orbital Zeeman 
term cx B z <j z t z couples B z to the pseudospin associated 
with the sublattice degree of freedom a z . However, the 
field B z does not induce a global pseudo-spin (sublattice) 
polarization, because this term has opposite signs in the 
two valleys so that the pseudospin polarization in the 
two valleys is antiparallel. As a result, the two sublat- 
tices remain indistinguishable even for finite B z . Only an 
electric field £ z polarizes the sublattices in BLG via the 
term cx £ z <t z tq, consistent with the fact that the pseu- 
dospin a z is even under time reversal.^ This behaviour 
is opposite to that of real spin, where B z s z Tq induces a 
real-spin polarization (s z ) (when averaging over the oc- 
cupied states in both valleys) while the term £ z s z t z does 
not, consistent with time reversal symmetry. Similarly, 
an in-plane magnetic field By gives rise to a macroscopic 
in-plane polarization of real spins, whereas the real-spin 
polarization induced by an electric field £n is anti-parallel 
in the two valleys. 



The fields By and £ y also couple to the in-plane pseu- 
dospin try. Term (6) of Table |l] induces an out-of-plane 
tilt of the spin orientation of individual states which, in 
the K' valley, has opposite signs for By and £y . Remark- 
ably, on average this yields an in-plane polarization (cry) 
which is nonetheless the same in each valley for fields 
By and £ y (Fig. [2j). This result reflects the fact that 
the macroscopic pseudospin polarization (cry) is neither 
even nor odd under time reversal.^ More precisely, in 
each valley the direction of (cry) is well-defined only up 
to a gauge-dependent angular offset. Yet the change of 
(cry) induced by a change in the in-plane orientation of 
the applied field is well-defined and it points clockwise in 
one valley and counterclockwise in the other valley (for 
both By and 5y and all terms in Tab. IT] giving rise to 
an in-plane pseudospin orientation of individual states). 
Specifically for the term (6), if we change the in-plane 
orientation of the fields By or £» by an angle ip, this 
changes the resulting average polarization (cry) by ±2cp. 
This implies, in particular, that reverting the direction 
of the external field yields the same orientation of the 
induced pseudospin polarization. We see here that the 
pseudospin polarization induced by external fields £ and 
B behaves qualitatively different from the polarization of 
real spins. 



Numeric values for the prefactors of the field-depen- 
dent interactions in Fig. [2] cannot be determined with 
the methods used here, but the important qualitative 
trends do not depend on these values. Therefore, in the 
calculations yielding Fig. [2] the prefactors were chosen 
large enough to provide a visually clear picture. 



B. Valley-contrasting axion electrodynamics 

Up to now, we have considered implications of the 
magneto-electric equivalence for interactions that are lin- 
ear in either B or £ . In every material, we also have in- 
teractions proportional to the squared field components. 
In SLG, these are interactions proportional to -By , B^, £,| 
and £ z . The magneto-electric equivalence implies that we 
may replace one factor in the square of one field by the 
corresponding components of the other field. We thus 
get scalar magneto-electric couplings 

c\\£\\-B\\t z + c z £ z B z t z , (3) 

as well as their trigonally anisotropic counterpart 

c A [{£ y B y - £ x B x )a x t z + (£ y B x + £ x B y )a y r ] . (4) 

Such scalars enter not only the Hamiltonian describing 
the dynamics of electrons in graphene, but they also en- 
ter the Lagrangian 18 from which we obtain Maxwell's 
equations in graphene. The bi-linear coupling of electric 
and magnetic fields displayed in Eq. Q is reminiscent 
of the £ ■ B contribution to the Lagrangian of axionsP 
hypothetical particles that were introduced to solve the 
strong CP problem in particle physics. Here the bi-linear 
couplings also have a valley dependence. Assuming in- 
dependent dynamics for the electrons in the two valleys, 
the new term in the Lagrangian becomes 

4x= J2 {-^)(c\\£wB {l +c z £ z B z ), (5) 

a=K,K' 

where = +1 and = — 1. The presence of £ ax 
modifies Gauss' law in BLG by inducing valley-dependent 
charge densities 

p {a) = v a [By • Vycy + (cy - C z ) Vy • By] (6) 

proportional to By, and it modifies Ampere's law by in- 
ducing current densities 

j( a) = v a z x [£ z V\\c z + (c z - cy) VyfJ (7) 

proportional to £ z (see Appendix|B]for details), thus real- 
izing a valley-contrasting axion electrodynamics in BLG, 
where the usual coupling 18 of charges to electric fields 
(Gauss' law) and currents to magnetic fields (Ampere's 
law) is completely reversed. 

The finiteness of any real sample leads to a spatial vari- 
ation of the coefficients Cj at its edge. As a result, equi- 
librium valley-isospin densities are induced at the BLG 
sheet's boundaries in the xy plane when the magnetic 
field has an in-plane component directed perpendicularly 
to such a boundary; and a perpendicular electric field 
generates equilibrium valley-isospin currents flowing par- 
allel to the system's boundaries in the xy plane. Figure[3] 
illustrates these effects in panels (a) and (b). The valley- 
asymmetric character of the edge densities and currents 
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E z <0 



fundamental research and new nano-spintronic applica- 
tions. Finally, it is fascinating that the phenomena dis- 
cussed here are present not only in bilayer graphene, but 
they exist in any material with similar symmetries. In 
particular, analogous behaviour c an be realized in ap- 
propriately designed metamaterials^^H or co fo atoms in 
optical lattices!^ The intrinsic tuneability of system pa- 
rameters in these latter realizations opens intriguing av- 
enues to tailor the relative magnitude of these matter- 
field interactions beyond the regime accessible in BLG. 



FIG. 3: Ramifications of the valley-contrasting axion electro- 
dynamics, (a) An in-plane magnetic field By induces equilib- 
rium valley-isospin densities p' Q ' at sample boundaries that 
are not aligned with B||. (b) A perpendicular electric field 

£ z gives rise to equilibrium valley-isospin edge currents j[i a '- 

(c) Valley-isospin densities can be created inside the sample 
by engineering inhomogeneity in the in-plane magnetic field. 

(d) Valley-isospin currents can be created inside the sample by 
engineering inhomogeneity in the perpendicular electric field. 



in BLG distinguishes them from similar edge effects in 
topological materials^ that generally require the pres- 
ence of strong spin-or bit c oupling. Also in contrast to 
the usually considered^ 21 cases where c z = c\\, valley- 
isospin densities (currents) can be induced in the bulk of 
BLG by spatially inhomogeneous fields By (£ z ). Scenar- 
ios for this are illustrated in panels (c) and (d) of Fig. |3j 
which can be realized, e.g., by placing nanomagnets on 
the surface of a BLG sample [for (c)] or arranging pairs 
of front and back gates with opposite voltage polarity 
around adjacent parts of the sample 2 ^ [for (d)]. More 
generally, the valley-contrasting axion electrodynamics 
discovered here enlarges the scope o f valle y-dependent 
electronic effects in graphene material d 21 * 22 ! by providing 
a comprehensive framework to design magneto-electric 
effects at boundaries and interface s in BLG, including 
associated valley-helical edge statesP^^ 
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Appendix A: Symmetry analysis and invariant 
expansion for BLG bandstructure: General 
discussion and results 

The Hamiltonian (JlJ describes the band structure of 
graphene as a function of the wave vecto r k m easured 
from the K point. The theory of invariant d 9 ! 14 ! makes it 
possible to also include perturbations 7C that are combi- 
nations of various quantities in addition to the wave vec- 
tor k, e.g., electric and magnetic fields S and B, strain 
e and the intrinsic spin s. The electron states at the K 
point transform according to the two-dimensional irre- 
ducible representation (IR) T 3 of the group D 3 (Ref HTJ) . 
The effective 2x2 Hamiltonian for the r 3 subspace can 
be expressed as 



III. DISCUSSION AND OUTLOOK 

One of the most exciting aspects of magneto-electric 
equivalence in BLG is the wide range of novel matter- 
field interactions resulting in this material. To achieve 
any type of effect within a valley, either an electric or a 
magnetic field component can be utilized. This includes, 
in particular, the different real-spin and pseudo-spin cou- 
plings of both Zeeman and spin-orbit character. Inge- 
nious combinations of nano-magnetic and nano-electronic 
fabrication capabilities can thus be used to realize new 
device architectures that are impossible to achieve in 
other materials. Additional variability is provided by 
the valley-isospin degree of freedom, allowing the selec- 
tive creation of global or valley-contrasting real-spin or 
pseudo-spin polarizations and their exploitation for both 



(Al) 



2 = 1 



Here a KX are prefactors, are 2x2 matrices that 

transform according to the IRs T K (of dimension L K ) con- 
tained in the product representation r 3 x r 3 = Ti +T 2 + 
r 3 of D3. Likewise, K, can be decomposed into irreducible 
tensor operators KS R that transform according to the 
IRs T K i of D 3 . Using the coordinate system in Fig. [I] we 
obtain the basis matrices and tensor operators listed in 
Tables HJ and EH 



Additional constraints for the Hamiltonian (Al) are 



due to time reversal invariance. The point group D$d of 
BLG contains symmetry elements R mapping the basis 
functions vI/ka at K on ^k'a a t K'. These basis functions 
are also mapped onto each other by the time-reversal 
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operation 9, i.e., we have 



(A2) 



with a u nitary matrix T , Combining these operations, 
we obtairp™ 



(A3) 



Here t denotes transposition and £ depends on the be- 
havior of K. under time reversal, k, B, and s are odd 
under time reversal so that then £ = — 1, while £ and e 
have £ = +1. Equation (A3) provides a general criterion 



for determining which terms in the expansion (All are 
allowed by time-reversal invariance and which terms are 
forbidden. The matrix T depends on the choice for the 
operation R. If R is the reflection R y at the yz plane 
[thus mapping the atoms in each sublattice in each layer 
onto each other, see Fig. [TJa)], we obtain 

(A4) 



H^{K) = H K {R~ 1 K) 



and the matrix T is simply the identity matrix. 

We note that under R y polar (p) and axial (a) vectors 
transform as 



Px -> -p x 



Py,z ^ Py,z 



(A5a) 
(A5b) 



The transformational properties for the components of 
the second-rank strain tensor can be expressed sim- 
ilarly. 14 If we denote by To (r z ) the unity (diagonal 
Pauli) matrix acting in valley- isospin space, Eq. (A5) 



can be summarized by writing general vector operators 
as (p x T z ,p y r ,p z T ) and (a x t , a y r z ,a z t z ). 

The group D3 characterizing the K point in BLG is a 
subgroup of the group D 3 h for the K point in SLGP^l so 
that any term allowed by spatial symmetries in "Hk(TC) 
for SLG is likewise allowed in BLG. Moreover, the con- 



straint ( A3 ) due to time reversal invariance is exactly 
equivalent to the constraint in SLGP Thus it follows im- 
mediately that the invariant expansion for BLG contains 
all terms that exist already for SLG. 

A more detailed analysis shows that the point group 
Dm for SLG distinguishes, as is usual, between polar vec- 
tors (such as the electric field £) and axial vectors (such 



TABLE II: Symmetrized matrices for the invariant expansion 
of the block for the point group D3. 

Block Representations Symmetrized matrices 



Hn 


ih x n = 


Ti 


Ih 


(i) 


H22 


r 2 xr; = 




Ih 


(i) 




iixr; = 


r 3 


r 3 


(M).(-M) 


-H23 


r 2 xr^ 


r 3 


r 3 


(1, -1), (-», -0 


-H33 


r 3 xr| 




Vi 


1 




= Li + r 2 


+ r 3 


r 2 
r 3 





as the magnetic field B). Thus each term in the SLG 
Hamiltonian with a certain functional form and linear in 
the field S or B is forbidden for the other field. However, 
the point group D3 c -D3;, relevant for BLG contains only 
rotations as symmetry elements so that it cannot distin- 
guish between polar and axial vectors. Therefore, the x 
and y components of any vector transform according to 
the IR F3 , whereas the z component transforms according 
to F2. This implies that spatial symmetries cannot dis- 
tinguish electric and magnetic fields in BLG. Moreover, 
Eq. (A3) treats electric and magnetic fields symmetri- 



cally, too. Thus it follows that every ^-dependent term 



TABLE III: Irreducible tensor components for the point 
group D3 (the group of the K point in BLG). Terms printed in 
bold give rise to invariants in Hk(K.) allowed by time-reversal 
invariance. (No terms proportional to k z are listed as they are 
irrelevant for graphene.) Contributions that are new in BLG 
(i.e., terms not allowed for the group D^h of the K point in 
SLG^are shown in red. Notation: {A, B} = \ (AB + BA). 

Fx I5 k x ~\~ ky j {Aia; j Sky k x } 5 B x k x ~t~ Byky ; 

+ £ + 

(^yy £xx ) k x -\- %€. X yky J Gyzk x € a 



-XZ ky 5 



€xx) B X ~\~ ^txyBy^ €y Z B X txZ By, 



{.^yy ^xx^&x ~(" 2i£xy£y 3 £yz&x ^xz&yi S x k x ~\~ Syky, 

&xB X -J- SyBy] S Z B Z \ S X Sx + 

i^S X ky Syk X ^)£zi SziJ^X^y ky£-x^)l 

Sxfeyy t-xx) ~\~ ^Sy€ X y j S x €y Z Sy€ xz 

\ky , 3k x ky } , B z , k x By kyB x ] £z i k x £ y ky£- x j 

SxBy — SyB X -) (^€ XX €yy^)ky -\~ *2,€ X yk X , Ey Z ky -\~ t XZ kxf 

i^-xx £yy} By ~\~ *2i€-xy B x j ( (~- x :v. ~\~ €yy) B z 3 tzzB z J 
£xzB x ~\~ tyzBy , [&xx £yy}£ y H - ^^-xy£ x 1 
( £xx ~\~ £yy^)£zi €zz£"Zi ^-xz^x 

&xky Sykxj S X By SyB X j S X £y SyEsXj 

{^S x k x ~\~ Syky^£ z , Sy (^€. X x ^yy) ~\~ "^^x^xy 1 
Sz^C-xx ~\~ Sx^xz Sytyz', Sz^-zz'i 

kx ) ky j { ky ~\- k X ) ky k X } ) 2{/Ca; , ky } , 

\kxi k x -\- ky\^ {fey? k x -\- ky\] 

B x , By , By ky B x k x , B x ky ~\- By k x , B z ky , B z k x , 

£x 7 £y 1 £y ky £/ X k x 1 &X ky ~\~ Ely kx 5 E ' Z ky , £■ Z kx 5 

£-y By £>x Bx ? £y B x ~\~ £.>' By , £- y B z ) £ X B z 5 



£■ z By , £■ z B x 



€-xx j *2i€-xy j Cyz ? 



( €-XX ~\~ ^yy) ( k -j. ? ky ) j €.yzk x ~\~ txzky ; €-xZ k x ^ 
( €.yy^)k x ~\~ %€. X yky , (€yy ^XX ) ky ~\~ "2C X yk x 

€-Zzkx) €-Zzky , 

(flS ~\~ £yy ) 5 By ) J £yzB X ~\- € XZ By , € X zB X €y Z By , 

(fra £yy S )B x -\- *2(Lxy By , {tyy t-xx) By -f" Qt€ X y B x , 
2£xy B z , (fxi €-yy^)B Zl € zz B Xj € zz By 

; £xzB z ,e yz B z ; 

{^xx H~ ^i/y ) (^a: 5 £y ) 5 £yz&x ~f~ ^-xz&y 1 €xz&x £yz&y 1 
(^xa; — tyy^j^x ~(" ^^xyEy 1 i_^yy ^.xx^&y H~ ^^-xy&x 1 
^^xy^z 1 (fsi )£z j £zz&x j ^zz£y 1 €-xz&z 1 €-yz&z 1 

Syky s x k x ^ s x ky -\- Sykx^ s z ky , $zkxi 

SyBy S X B X }S X By + SyB 

X ; S zBy ; S Z B X 1 

s„B z ,- - S X B z ; Sy£y Sx^x) S x £y + 
Sz£y<> Sz^x'i Sy&zi &x£ z] 

&z (kx£ y ~\~ ky£>x ) j ( k x £- x kyS- y} 1 
Syk x }£z) (^Sxkx Syky 

(Sx; H~ ^-yy^) 1 ^Sz^-xy •> Szi^-xx ^yy^)i 

Sxi^-XX ~~ ^yy^) ^Sy€ X y , Sy(^€yy fra) 2Sa;£a;yi 

5a;£zZj Sy^-ZZl Sz^-XZ) S Z €y Z ^ S X €y Z -f - Syt XZ -> S X £ XZ SyEy; 
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in the BLG Hamiltonian ( Al I is accompanied by another 



term where £ is simply replaced by B (and vice versa for 
B-dependent terms). However, the prefactors of these 
terms are, in general, unrelated .™ Table [i] summarizes 
some of the new terms arising from this magneto-electric 
equivalence. 

Similar to SLG, pseudospin up and down in BLG cor- 
responds to atoms in sublattice B and B', see Fig. [ija) .™ 
With the convention that for both valleys v = K, K' the 
up and down eigenstates of pseudospin er 2 correspond to 
the same sublattice, the above procedure determines the 
Hamiltonians rl v up to phases <j) v corresponding to pseu- 
dospin rotations <jj = exp(i(/) v a z /2) <jj exp(— i(j) v a z /2) of 
the basis matrices Oj (j = x, y) entering the invariant ex- 



pansion (Al) (Ref I50|) . Due to this U(l) gauge freedom 



the direction of the average in-plane spin polarization 
(<Ti|)„ induced within a valley by fields B and £ is well- 
defined only up to a uniform field-independent angular 
offset 6„. 



Appendix B: Valley-contrasting axion 
electrodynamics in BLG: Calculation of induced 
densities and currents 

The presence of £ ax leads to modifications of the 
inhomogeneous Maxwell's equations (i.e., Gauss' and 
Ampere's laws) amounting to the replacements^ 

D ->■ D - (c|| B|| + c z B z z) , (Bla) 

a 

H -> H + ^i/ a ( C || 5|| +c z £ z z) , (Bib) 

a 

where D and H denote the macrosco pic e lectric and mag- 
netic fields in matter .™ Equations (Bl| imply the exis- 



tence of extra valley-dependent charge densities and cur- 
rents that are most generally given by 



P {a) = v*V- (c|,B|,+c„B,S), 
j (Q) = -v a V x (cy£|| +c z £ z z) 
-v a d t (en By +c z B z z) . 



(B2a) 



(B2b) 



The coefficients Cj arc finite constants within BLG but 
must vanish outside the sample. Hence, there is a bound- 
ary contribution to Vc, = VyCj, and we have 



' ex [(V,| C ||) • B,| + C || V • B,| + c z d z B z ] . (B3) 



Using the homogeneous Maxwell equation V • B = 0, the 
magnetic-field-induced charge densities are then given by 



Eq. Furthermore, the two-dimensionality of BLG 

implies that currents flow only in-plane. Hence, both 
(V||C||) x f || and d t B z z do not contribute to the current, 
which is just given by 



■ la) 



[ C || (Vx£||+d t B||)+Vx(c 2 £ 2 z)] 



(B4) 

The other homogeneous Maxwell equation can be written 
as V x f || + <9fB|| = —V x (£ z z) — dt B z z whose in- 
plane projection inserted into Eq. (B4) yields Eq. (JtJ) . 



For homogeneous, time-independent fields £ and B the 
ordinary couplings of £ and B to p and j vanish so that 
we only get the valley-contrasting axion electrodynamics 
of Eqs. @ and 0. 

The valley-contrasting axion electrodynamics dis- 
cussed here can be expected to apply in situations when 
the valley-isospin degree of freedom is conserved. This 
will generally be the case for inhomogeneities induced by 
smoothly varying external fields (such as the one pro- 
posed in Ref l2H|) . In contrast, the physical termination 
of the BLG sheet is abrupt on an atomic scale, and de- 
tails will usually matter™ Nevertheless, the recently ob- 
served™ emergence of a universal behavior at disordered 
BLG edges consistent with that obtained from our valley- 
contrasting axion-electrodynamics suggests a wider ap- 
plicability of this formalism. 



Appendix C: Magnitude of Prefactors in Eq. ([2]) 

Recent experiments™ demonstrated that a displace- 
ment field of ~ 1 V/nm generates a bandgap of ~ 0.1 eV 
in BLG, implying that X z ~ 0.1 enm (with electron 
charge e) consistent with first-principles calculations.™ 
For the remaining prefactors in Eq. ([2]), the Slonczewski- 
Weiss-McClure (SWM) modeP applied to BLG yields 
(see also Ref lTTj) 



47074 



6.2 



7l7a 

7o + 7! Mb A 2 



7? 



2 7a 



(CI) 



3xlO~ 4 enm/T. (C2) 



The SWM parameters 70, 71, and 74 are defined in 
Refl2l)l and we have -f a — 2h 2 /(3m a 2 ), where mo is the 
electron mass in vacuum and a is the lattice constant of 
BLG's planar honeycomb structure. 
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